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0^ Abstract 

(N 

jJere we consider the possibility— envisaged by many authors as feasible in the near future— of measuring at 10% or better the moment 
I iqf inertia / of the pulsar J0737-3039A via the gravitomagnetic spin-orbit periastron precession (analogous to the Lense-Thirring 
Jn^^ericentre precession in the case of a test-particle orbiting a central spinning mass) . Such a gravitomagnetic effect is expected to be 
I of the order of lO"'' deg yr~^ and the present-day precision in measuring the periastron precession of J0737-3039A via pulsar timing 
^ is 6.8 X 10^* deg yr~^. However the systematic uncertainty in the much larger first-order post-Newtonian (IPN) gravitoelectric 
i_^precession (analogous to the Einstein Mercury's perihelion precession in the weak- field and slow-motion approximation), which 
should be subtracted from the measured one in order to pick up the gravitomagnetic rate, is of primary importance. Indeed, by 
00 determining the sum of the masses by means of the third Kepler law, such a bias amounts to 0.03165 deg yr^^, according to the 
^ fcurrent level of accuracy in knowing the parameters of the J0737-3039 system. The major sources of uncertainty are the Keplerian 
r~ projected semimajor axis xb of the component B and the post-Keplerian parameter s, identified with sini; their knowledge should 
be improved by three orders of magnitude at least; the bias due to the Keplerian projected semimajor axis xa of the component 
A amounts to ~ 10% today. The present-day level of accuracy in the eccentricity e would affect the investigated measurement at a 
' \ percent level, while the impact of the orbital period Pb is completely negligible. If, instead, the sum of the masses is measured by 
lO means of the post-Keplerian parameters r and s, it turns out that r should be measured five orders of magnitude better than now: 
according to the present-day level of accuracy, the total uncertainty in the IPN periastron rate is, in this case, 2.11819 deg yr^^. In 
conclusion, the prospect of measuring the moment of inertia of PSR J0737-3039A at 10% accuracy or better seems unlikely given the 
limitations to the precision with which the system's basic binary and post-Keplerian parameters can be measured via radio timing. 
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1. Introduction 

The measurement of the moment of inertia / of a neutron 
star at a 10% level of accuracy or better would be of crucial 
importance for effectively constraining the Equation-Of- 
Statc (EOS) describing matter inside neutron stars (Mor- 
rison et al., 2004; Bejger et al., 2005; Lattinicr and Schutz, 
2005; Lavagetto et al, 2007). 

After the discovery of the double pulsar PSR J0737- 
3039A/B system (Burgay et al., 2003), whose relevant or- 
bital parameters are listed in Table 1, it was often argued 
that such a measurement for the A pulsar via the post- 
Newtonian gravitomagnetic spin-orbit periastron preces- 
sion (Barker and O'Connell, 1975a; Damour and Schaefer, 
1988; Wex, 1995) would be possible after some years of 
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accurate and continuous timing. Lyne et al. (2004) write: 
"Deviations from the value predicted by general relativity 
may be caused by contributions from spin-orbit coupling 
(Barker and O'Connell, 1975b), which is about an order of 
magnitude larger than for PSR B1913-I-16. This potentially 
will allow us to measure the moment of inertia of a neutron 
star for the first time (Damour and Schaefer, 1988; Wex, 
1995)." 

According to Lattimer and Schutz (2005), "measurement 
of the spin-orbit perihelion advance seems possible." 

In (Kramer et al., 2006) we find: "A future determination 
of the system geometry and the measurement of two other 
PK parameters at a level of precision similar to that for 
w, would allow us to measure the moment of inertia of 
a neutron star for the first time (Damour and Schaefer, 
1988; Wex, 1995). [...] this measurement is potentially very 
difficult [...] The double pulsar [...] would also give insight 
into the nature of super-dense matter." 



In (Damour, 2007) it is written: "It was then pointed out 
(Damour and Schacfer, 1988) that this gives, in principle, 
and indirect way of measuring the moment of inertia of 
neutron stars [...]. However, this can be done only if one 
measures, besides ^ k, two other PK parameters with 10^^ 
accuracy. A rather tall order which will be a challenge to 
meet." 

Some more details are released by Kramer et al. (2005): 
"[...] a potential measurement of this effect allows the mo- 
ment of inertia of a neutron star to be determined for the 
first time (Damour and Schacfer, 1988). If two parameters, 
e.g. the Shapiro parameter s and the mass ratio i?, can 
be measured sufficiently accurate, an expected t^cxp can be 
computed from the intersection point." 

Here we wish to examine, more precisely and with more 
quantitative details than done in the existing literature, the 
conditions which would make feasible to measure /a at 10% 
or better in the PSR J0737-3039A/B system in view of the 
latest timing results (Kramer et al., 2006). In particular, we 
will show how crucial the impact of the mismodelling in the 
known precessional effects affecting the periastron rate of 
PSR J0737-3039A/B is if other effects are to be extracted 
from such a post-Keplerian parameter. Such an analysis will 
turn out to be useful also for purposes other than measuring 
gravitomagnetism like, e.g., putting more severe and real- 
istic constraints on the parameters entering various mod- 
els of modified gravity. Indeed, in doing that for, e.g., a 
uniform cosmological constant A Jetzer and Sereno (2006) 
took into account only the least-square covariance sigma 
of the estimated periastron rate (6.8 x lO^"' deg yr~^): the 
systematic bias due to the first post-Newtonian (IPN) pe- 
riastron precession was neglected. 

2. The systematic uncertainty in the IPN 
periastron precession 

By assuming / sa 10^^ kg m^ (Morrison et al., 2004; Be- 
jger et al., 2005), the gravitomagnetic spin-orbit periastron 
precession is about ojgm ~ 10~^ deg yr~^, while the error 
Jwnioas with which the periastron rate is phenomcnologi- 
cally estimated from timing data is currently 6.8 x 10"'* 
deg yr~^ (Kramer et al., 2006). In order to measure the 
gravitomagnetic effect— and, in principle, any other dynam- 
ical feature affecting the periastron— <^Wmcas is certainly 
of primary importance, but it is not the only source of 
error to be carefully considered: indeed, there are other 
terms contributing to the periastron precession (first and 
second post-Newtonian, quadrupole, spin-spin (Barker and 
O'ConneU, 1975a; Damour and Schacfer, 1988; Wex, 1995)) 
which must be subtracted from Wmcas, thus introducing a 
further systematic uncertainty due to the propagation of 
the errors in the system's parameters entering their analyt- 
ical expressions. A preliminary analysis of such aspects, can 
be foimd in (Lattimer and Schutz, 2005). However, apart 
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from the fact that its authors make use of the value for i 
measured with the scintillation method (Coles et al., 2005) 
which is highly uncertain for the reasons given below, in 
using the third Kepler law to determine the sum of the 
masses they also confound the relative projected semima- 
jor axis^ asini (see eq. (3)) with the barycentric projected 
semimajor axis x, which is the true measurable quantity 
from timing data, so that their analysis cannot be consid- 
ered reliable. The semimajor axis a of the relative motion 
of A with respect to B in a binary system is just the sum 
of the semimajor axes Cbc of A and B with respect to the 
system's barycenter, i.e. a = a^^ + a^^. 

Let us, now, consider the largest contribution to the peri- 
astron rate, i.e. the IPN precession (Damour and DerucUe, 
1986; Damour and Taylor, 1992) 

where Tq = GMq / and M = niA + , in Solar masses, . 
It is often referred to as gravitoelectric in the weak-field 
and slow-motion approximation: in the context of the Solar 
System it is the well known Einstein Mercury's perihelion 
precession of about 43 arcsec cy~^. Thus, 

(2) 

The sum of the masses M enters eq. (1); as we will see, 
this implies that the relative semimajor axis a is required 
as well. For consistency reasons, the values of such param- 
eters used to calculate eq. (1) should have been obtained 
independently of the periastron rate itself. We will show 
that, in the case of PSR J0737-3039A/B, it is possible. 
Let us start from the relative semimajor axis 

a={l + R]( ^) = 8.78949386 x 10® m. (3) 
Vsmi/ 

It is built in terms of R, the projected semimajor axis xa 
and sin i; the phenomenologically estimated post-Keplerian 
parameter s determining the shape of the logarithmic 
Shapiro time delay can be identified with sini in general 
relativity and the ratio R = xb/ xa has been obtained from 
the phenomenologically determined projected semimajor 
axes, being equal to the ratio of the masses in any Lorentz- 
invariant theory of gravity (Damour and DerucUe, 1985; 
Damour and Schacfer, 1988; Damour and Taylor, 1992) 

The uncertainty in a can be conservatively evaluated as 
5a < 6a\R + Sa\s + Sa\xj^, (5) 

^ Our a must not be confused with the one used in (Lattimer and 
Schutz, 2005) which is, in fact, a^c- 



with 

Sa\R < (—) SR = 466, 758 m, 
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Sals < a I — ) = 342,879 m. 



(6) 



5a\ 
Thus 



..<a(^^)=621m. 



Sa < 810,259 m. 

eq. (7) yields a relative uncertainty of 

5a . 
— = 9 X 10"^. 
a 



(7) 



(8) 



It is important to note that zb, via i?, and s have a major 
impact on the overall uncertainty in a; our estimate has to 
be considered as conservative because we adopted for bs 
the largest value quoted in (Kramer et al., 2006). In regard 
to the inclination, we did not use the more precise value for 
i obtained from scintillation measurements in ^ (Coles et 
al., 2005) because it is inconsistent with that derived from 
timing measurements (Kramer et al., 2006). Moreover, the 
scintillation method is model-dependent and it is not only 
based on a number of assumptions about the interstellar 
medium, but it is also much more easily affected by various 
other effects. However, we will see that also xa has a non- 
negligible impact. Finally, let us note that we purposely 
linearly summed up the individual sources of errors in view 
of the existing correlations among the various estimated 
parameters (Kramer et al., 2006). 

Let us, now, determine the sum of the masses: recall that 
it must not come from the periastron rate itself. One possi- 
bility is to use the phenomenologically determined orbital 
period Pb and the third Kepler law getting ^ 



GM = — 



27r 



(9) 



With eq. (3) and eq. (9) we can, now, consistently calculate 
eq. (1) getting 



WlPN 



(l-e2) 



XA + Xb 



27rY 



16.90410 deg yr" 



in this way the IPN periastron precession is written in 
terms of the four Keplerian parameters Pb, e, xa, ^cb and of 
the post-Keplerian parameter s. The mismodeling in them 
yields 



^ See also (Lyutikov and Thompson , 2005). 

^ In principle, also the IPN correction to the third Kepler law cal- 
culated by (Damour and Deruelle, 1986) should be included, but it 
does not change the error estimate presented here. 



<5wipnUb — 2cjiPN 



5x\ 



(a;A + xb) 



0.01845 degyr-\ 



(^wipnIs < 2[jiPN ( — ) = 0.01318 deg yr \ 



< (5wipnUa ^ SwipN 



^wipnU < 2ea;iPN 



(xA -f a;B) 



5e 
(1^ 



= 1 X 10"^ degyr-\ (11) 



2 X 10"^ deg yr"\ 



(5wiPN|Pb < St^iPN ( ) = 0(10 ^) deg yr \ 



Thus, the total uncertainty is 
(SwiPN < 0.03165 deg yr"\ (12) 
which maps into a relative uncertainty of 

^ = 1.8 X 10-=^. (13) 

O^lPN 

As a consequence, we have the important result 

Atj = tj^eas - tl'iPN = (-0.00463 ± 0.03233) deg yr"^(14) 

Every attempt to measure or constrain effects predicted by 
known Newtonian and post-Newtonian physics (like, e.g., 
the action of the quadrupole mass moment or the gravit- 
omagnetic field), or by modified models of gravity, for the 
periastron of the PSR J0737-3039A/B system must face 
with the bound of eq. (14). 

Should we decide to use both the post-Keplerian param- 
eters related to the Shapiro delay (Damour and Deruelle, 
1986; Damour and Taylor, 1992) 



2tt 



-2/3 



(15) 



for determining the sum of the masses, we would have, with 



WiPN — 



n X 3/2 / X 9/4 
Pb\ I r ^ ' 



,19/4 



(l-e2)V27ry \xa) {xA + x^f/'^ 
which yields 

cOiPN = 17.25122 ± 2.11819 deg yr'^ 



(16) 



(17) 



The major source of uncertainty is r, with 2.06264 deg yr"'^; 
the bias due to the other parameters is about the same as 
in the previous case. 



Q 



Let us, now, consider the second post-Newtonian contri- 
bution to the pcriastron precession (Damour and Schaefer, 
1988; Wex, 1995) 



(jJ2PN 



3(GAf)^/2 



c4aV2(l 



13 / m\ + mg\ 32 (mpjn^ 



AP 



3 V AP 



up to terms of order 0{e^). For our system it amounts to 
4 X lO^'* deg yr^^, so that it should be taken into account 
in Aiu. However, it can be shown that the bias induced 
by the errors in AI and a amounts to 4 x 10~^ deg yr~^, 
thus affecting the gravitomagnetic precession at the percent 
level. 



skepticism than done so far. Measuring gravitomagnetism 
is a challenging enterprizc. 
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3. Discussion and conclusions 

O'Connell (2004), aware of the presence of other non- 
gravitomagnetic contributions to the pulsar's pcriastron 
rate, proposed to try to measure the gravitomagnetic 
spin-orbit precession of the orbital angular momentum 
(Barker and O'Connell, 1975a) (analogous to the Lense- 
Thirring node precession in the limit of a test particle 
orbiting a massive body) which is not affected by larger 
gravitoelectric contributions. However, its magnitude is 
~ (10~^ deg yr"^) sinip^ where ip is the angle between the 
orbital angular momentum and the pulsar's spin; thus, it 
would be negligible in the PSR J0737-3039A/B system be- 
cause of the near alignment between such vectors (Stairs 
et al., 2006), in agreement with the observed lack of profile 
variations (Manchester et al., 2005; Kramer et al., 2006). 

In regard to the measurement of the moment of iner- 
tia of the component A via the gravitomagnetic pcrias- 
tron precession, our analysis has pointed out that the bias 
due to the mismodelling in the IPN gravitoelectric contri- 
bution to pcriastron precession-expressed in terms of the 
phenomenologically measured parameters Ph-, e, xa, xb-, s- 
is the most important systematic error exceeding the ex- 
pected gravitomagnetic rate, at present, by two orders of 
magnitude; the major sources of uncertainty in it are xb 
and s, which should be measured three orders of magni- 
tude better than now to reach the 10% goal. The projected 
semimajor axis xa of A, if known one order of magnitude 
better than now, would induce a percent-level bias. Instead, 
expressing the IPN gravitoelectric pcriastron rate in terms 
of Pb, e, xa, a;B, s, f would be definitely not competitive be- 
cause the improvement required for r would amount to five 
orders of magnitude at least. We prefer not to speculate 
now about the size of the improvements in timing of the 
PSR J0737-3039A/B system which could be achieved in fu- 
ture. Since the timing data of B are required as well for xb 
and in view of the fact that B appears as a strong radio 
source only for two intervals, each of about 10-min dura- 
tion, while its pulsed emission is rather weak or even un- 
detectable for most of the remainder of the orbit (Lyne et 
al., 2004; Burgay et al., 2005), the possibility of reaching in 
a near future the required accuracy to effectively constrain 
/a to 10% level or better should be considered with more 
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Table 1 

Relevant Kcplcrian and post-Keplerian parameters of the binary sys- 
tem PSR J0737-3039A/B (Kramer et al., 2006). The orbital period 
Plj is measured with a precision of 4 X 10~® s. The projected semima- 
jor axis is defined a.s x = (ai,c/c) sini, where ate is the barycentric 
semimajor axis (the relative semimajor axis a = (xa + x^)c/ sini), 
c is the speed of light and i is the angle between the plane of the 
sky, perpendicular to the line-of-sight, and the orbital plane. The ec- 
centricity is e. The best determined post-Keplerian parameter is, to 
date, the periastron rate u! of the component A. The phenomenolog- 
ically determined post-Keplerian parameter s, related to the general 
relativistic Shapiro time delay, is equal to sini; we have conserva- 
tively quoted the largest error in s reported in (Kramer et al., 2006). 
The other post-Keplerian parameter related to the Shapiro delay, 
which is used in the text, is r. 

Pb (d) XA (s) XB (s) e uj (deg yr~^) s r (^ts) 

0.10225156248(5) 1.415032(1) 1.5161(16) 0.0877775(9) 16.89947(68) 0.99974(39) 6.21(33) 



